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Abstract 

First-order  nonlinear  interactions  of  Tollmien-Schlichting  waves  of 
different  frequencies  and  initial  amplitudes  in  boundary- layer  flows  are 
analyzed  by  using  the  method  of  multiple  scales.  For  the  case  of  two 
waves,  a strong  nonlinear  interaction  exists  if  one  of  the  frequencies  3 
u)2^s  twice  the  other  frequency.ldTT^Iumerical  results  for  flow  past  a 
flat  plate  show  that  this  interaction  mechanism  is  strongly  destabilizing 
even  in  regions  where  either  the  fundamental  or  its  harmonic  is  damped 
in  the  absence  of  the  interaction.^  For  the  case  of  three  waves,  a 
strong  nonlinear  interaction  exists  when  u)3  = u)2-  on.  This  combination 
resonance  causes  the  amplitude  of  the  wave  with  the  difference  frequency 
o)3  to  multiply  many  times  in  magnitude  in  a short  distance  even  if  it 
is  damped  in  the  absence  of  the  interaction.  The  initial  amplitudes 
play  a dominant  role  in  determining  the  changes  in  the  amplitudes  of 
the  waves  in  both  of  these  mechanisms. 


I.  INTRODUCTION 


* 


One  of  the  major  roads  from  laminar  to  turbulent  flow  involves  the 
initial  linear  amplification  of  disturbances,  which  might  be  present  in 
the  flow.  However,  as  these  disturbances  grow  to  appreciable  amplitudes, 
nonlinear  effects  set  in.  The  nonlinear  mechanisms  that  are  activated 
depend  on  the  spectrum  of  the  disturbances.  In  this  paper,  we  investigate 
two  of  these  mechanisms. 

In  his  experiments  on  the  transition  from  laminar  to  turbulent  flow 

in  a separated  shear  layer,  Sato1  observed  the  appearance  of  the  subharmonic 

of  order  one-half  in  addition  to  the  higher  harmonics  of  the  fundamental 
2 

wave.  Wille  observed  the  development  of  subharmonic  waves  while  investigating 

3 

the  stability  of  both  circular  and  plane  jets.  Kachanov,  et  al  observed 

that,  in  addition  to  the  higher  harmonics  of  a fundamental  wave,  which 

was  introduced  in  the  flow  by  a vibrating  ribbon,  a subharmonic  wave 

with  one-half  the  frequency  of  the  fundamental  wave  appeared  downstream. 

4 

Michalke  postulated  that  the  subharmonic  appears  when  two  vortices 

5 

rotate  around  each  other  in  a fusion  mating  dance.  Kelly  showed  that 
the  appearance  of  the  subharmonic  in  a shear  layer  is  due  to  a secondary 
linear  instability  associated  with  a time-dependent  flow  that  consists 
of  the  superposition  of  the  basic  flow  and  a finite-amplitude  funda- 
mental wave.  Nayfeh  and  Bozatli5 6  investigated  the  appearance  of  the 
subharmonic  in  boundary  layers  by  analyzing  the  instability  associated 
with  a time-dependent  flow  that  consists  of  the  superposition  of  the 
basic  flow  and  a Tollmien-Schlichting  wave.  The  results  show  that  the 
amplitude  of  the  fundamental  wave  must  exceed  a critical  value  to 
trigger  this  parametric  instability.  This  value  is  proportional  to  a 
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detuning  parameter  that  is  the  real  part  of  k - 2K,  where  k and  K are 

the  wavenumbers  of  the  fundamental  and  its  subharmonic,  respectively. 

For  the  Blasius  flow,  the  critical  amplitude  is  approximately  29%  of  the 

mean  flow.  For  other  flows  where  the  detuning  parameter  is  small,  such 

as  free-shear  layer  flows,  the  critical  amplitude  can  be  small,  thus  the 

parametric  instability  might  play  a greater  role.  Since  the  analysis  of 
5 6 

Kelly  and  Nayfeh  and  Bozatli  are  linear,  they  do  not  account  for  the 
effect  of  the  subharmonic  wave  on  the  fundamental  wave.  This  effect  may 
be  small  initially,  but  as  the  subharmonic  grows  appreciably,  its 
effect  on  the  fundamental  cannot  be  neglected.  One  of  the  purposes  of 
the  present  paper  is  to  determine  the  nonlinear  interaction  of  a Tollmien- 
Schlichting  wave  with  its  subharmonic. 

7 0 Q 

Sato  , Miksad  , and  Kachanov,  et  al  observed  that  the  nonlinear 
development  of  the  waves  in  the  transition  region  depends  on  the  initial 
and  external  disturbances.  Sato7  conducted  an  experiment  on  the  stability 
of  symmetric  laminar  wakes  by  exciting  two  unstable  modes  with  the 
frequencies  fi  and  f2.  He  observed  the  generation  of  waves  having  the 

O 

frequencies  f2  ± f i . Miksad  excited  two  unstable  modes  of  a laminar 

asymmetric  free-shear  layer.  He  also  observed  nonlinear  triggered 

instabilities  of  the  difference  mode  f2  - fi,  subharmonics,  and  higher 

g 

harmonics  of  the  fundamental  waves.  Kachanov,  et  al  introduced  two 
Tollmien-Schlichting  waves  in  the  boundary  layer  on  a flat  plate  b> 
using  two  vibrating  ribbons.  They  observed  the  appearance  and  growth  of 
a Tollmien-Schlichting  wave  having  the  difference  frequency  f2  - fi. 
Norman^  also  observed  the  amplification  of  the  difference  harmonic  of 
two  introduced  disturbance  waves  in  his  experimental  study  of  secondary 


3 


flows  around  and  downstream  of  protuberances  in  laminar  boundary  layers. 
The  second  purpose  of  the  present  paper  is  to  determine  the  nonlinear 
interaction  of  three  Tollmien-Schl ichting  waves  (combination  resonance) 
in  boundary  layers  and  show  that  the  difference  frequency  can  be  very 
unstable  when  generated  by  the  nonl inearity,  even  though  it  is  stable 
when  introduced  by  itself  in  the  boundary  layer. 

The  problem  is  formulated  in  Sec.  II.  The  analysis  for  the  com- 
bination-resonance case  is  contained  in  Sec.  Ill,  while  the  results  for 
the  second-harmonic  case  are  stated  in  Sec.  IV.  The  numerical  procedure 
is  discussed  in  Sec.  V,  while  the  numerical  results  are  presented  in 
Sec.  VI. 
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II.  PROBLEM  FORMULATION 


We  consider  nonlinear  interactions  of  wave  packets  in  a two-dimensional 
steady  incompressible  boundary -layer.  The  equations  describing  the 
motion  of  the  fluid  are 


sQ 

3x  3y  u* 

3a  + 0iy.+  v^  = -^  + iv2u 

3t  u 3x  v 3y  3x  IT  v u* 


22  + op  + V^  - - 22  + lv2v 
3t  “5x  3y  3y  R v v* 


u = v = 0 at  y = 0. 


(1) 

(2) 

(3) 

(4) 


u ■+  1 


as 


(5) 


where 


V2  = 


32 

Sx2 


W 


Here,  x and  y are  made  dimensionless  by  using  a reference  length  6r, 
the  time  is  made  dimensionless  by  using  <S  /U  , and  the  velocities  are 
made  dimensionless  by  using  the  freestream  velocity  U^.  The  Reynolds 
number  R = U^/v  with  v being  the  fluid  kinematic  viscosity. 

The  analysis  is  restricted  to  basic  flows  that  are  slightly 
nonparallel  ( i . e. , vary  slowly  in  the  streamwise  direction)  and  to 
disturbances  that  are  small  but  finite.  The  slow  variation  is  expressed 
by  using  the  slow  scale  xi  seix,  where  ei  is  a small  dimensionless  quantity 
that  characterizes  the  nonparallelism  of  the  flow  and  can  be  related 
to  R by  ei  = R";.  The  smallness  of  the  amplitude  of  the  disturbance 
is  expressed  by  introducing  the  small  dimensionless  parameter  e.  For 
a general  solution,  we  assume  that  e = 0( e x ) so  that  the  resulting 
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expansion  accounts  simultaneously  for  the  effects  of  nonparallelism  and 
nonl inearity.  When  e « ei,  the  nonlinear  effects  are  negligible  and 
the  solution  reduces  to  those  obtained  in  Refs.  11  and  12.  When 
e » ei,  the  nonparallel  effects  are  negligible  and  the  solution  reduces 
to  equations  with  constant  coefficients. 

We  assume  that  each  flow  quantity  is  the  sum  of  a mean-flow  quan- 
tity and  an  unsteady  disturbance  quantity,  which  is  assumed  to  be  much 
smaller  than  the  mean-flow  quantity.  We  can  then  express  the  velocity 
components  and  the  pressure  as 

u(x,y,t)  = U0(xi,y)  + eu(x,y,t) , (6) 

v(x,y,t)  = eiV0(xi,y)  + ev(x,y,t),  (7) 

p(x,y,t)  = P (xj  + ep(x,y,t),  (8) 


where  U0,  V0,  and  P0  are  the  nonparallel  basic-flow  quantities.  Sub- 
stituting Eqs.  (6)  - (8)  into  Eqs.  (1)  - (5)  and  subtracting  the  basic- 
flow  quantities,  we  obtain 

§7*1?  * °-  <9> 


3u  3u 
eU  3x  ' £V  37 


iv  + u 3V  3£  1 2 

3t  0 3x  3y  R 


2 3V0  ,,  3v  3Vr 

ElU  3Xj  " ElV  3y  " ElV  3y 


..  „ 3v 

EU  37  • £V  3^  ’ 


u = v = 0 at  y s 0, 


u,  v as  y 


1 


Without  loss  of  generality,  we  let  e = ej.  To  determine  the  wave- 

packet  solutions  of  Eqs.  (9)  - (13),  we  use  the  method  of  multiple 
13 

scales  and  seek  an  expansion  in  the  form 

u = Uo(x0,Xi,y,T0,Ti)  + eUi (x0 ,Xi ,y,T0 ,Ti ) + ....  (14) 

v = Vo(x0,Xi ,y,T0,Ti)  + eVi(x0,Xi,y,T0,Ti)  + ....  (15) 

P = Po(x0,Xi,y,T0,Ti)  + cpi(xo,Xi,y,To,Ti)  + ....  (16) 

where  x<j  = x,  T0  = t,  and  Ti  = et.  Substituting  Eqs.  (14)  - (16)  into 
Eqs.  (9)  - (13)  and  equating  coefficients  of  like  powers  of  e,  we  obtain 


Order  e° 


9uo  , 9v0  _ „ 

9x0  ay 

i (uo ,vo ,p0 ) = Irt  + Uo  ISt + Vo  ly9"  + ^ v°Uo 

2(uo,vo,Po)S  Irt  + 1,0  + !?MV°V<>  = 0’ 


u0  = v0  = 0 

at 

y = 0 

u o , Vo  -*■  0 

as 

y 00 

Order  e 

9Uj.  + 

9vi 

3uo 

9x0 

9y 

9Xi 

0, 


l(Ul,V!,Pi) 


9uo  _ .I  3up  _ 3Po  . 2 92u  p 
9T i ® 9Xi  9Xj  R"  9Xq3Xi 


- Vo 


9Up 

9y 


Uo 


9u  o 
9xo 


' ' ^2  (Ul.Vi  ,pl  ) = - 1^-  - 


■ v0 

9Vo 

9X) 


9u0 

9y 


+ 1 

R 3x09Xi 


92v0 


_ u ?y.ji  _ v 
V°  3xo  0 


9 )Li 
9y 


(17) 

(18) 

(19) 

(20) 
(21) 

(22) 

(23) 

(24) 
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In  what  follows,  we  describe  the  details  of  the  analysis  for  the 
combination-resonance  case  and  only  state  the  results  for  the  second- 
haromnic  resonance  case. 


III.  COMBINATION  RESONANCES 


t 


, 


A.  First-Order  Problem 

For  the  case  of  combination  resonances,  we  consider  three  wave- 
packets  centered  at  the  frequencies  on,  uj2,  and  u>3,  Then  we  examine  the 
resonances  that  might  exist  among  them.  Thus,  the  solution  of  Eqs. 

(17)  - (21)  is  expressed  as  a linear  combination  of  three  Tollmien- 
Schlichting  waves;  that  is,. 

u0=  A1(xi,T1)cu(y;xi)exp(i6i)  + A2(xi,Ti)  x 

Ci2(y;xi)exp(i02)  + A3(xi,Ti)d3(y;xi)exp(i03)  + c.c.,  (27) 
v0=  Ai(xi,Ti)c2i(y;xi)exp(i0i)  + A2(xi,Ti)  x 

C22(y;xi)exp(i02)  + A3(xi,Ti)c23(y;xi)exp(i03)  + c.c.,  (28) 

Po=  Ai(x1,T1)c31(y;Xi)exp(i0i)  + A2(x3,Ti)  x 

C32(y;xi )exp(i02)  + A3 (xi ,Ti )c33 (y;xi )exp(i03 ) + c.c.,  (29) 

where 

30„  80„ 

= kn(xi)’  aT7=-a)n  (n  = 1.  2.  3)  (30) 

with  the  being  real  constants.  The  quasi-parallel  Orr-Sommerfeld 
problems  for  these  waves  are 

kn>  s % + ik,c,"  ' °'  (3,) 

Sn>  c,n;  kn-  “n>  5 1(u»kn  * + V1* 

'32> 
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c2n’  c3n’  kn’  “Vi*  ~ 1(U°kn  ' wn^2n  + °c3- 


(°2  - k^2n  " 


C n = C n = 0 
in  2n 


y = o. 


V ^n’  0 as 


where  D = 3/3 y. 


B.  Second-Order  Problem 


Substituting  Eqs.  (27)  - (29)  into  Eqs.  (22)  - (26),  we  find  that 
the  inhomogeneous  parts  in  Eqs.  (22)  - (26)  contain  terms  proportional 


exp(i6i),  exp(i02),  exp(ie3), 

exp[i (02  - 0i)],  exp[i(0i  + 03)],  exp[i(02  - 03)] 

where  the  overbar  indicates  the  complex  conjugate.  The  terms  that  are 
proportional  to  these  exponential  expressions  will  create  secular  terms 
in  the  particular  solutions  for  ui,  v i,  and  pj  if  k3  ~ k2  - ki 
and  u3  ~ (d2  - oh;  that  is,  when  a combination  resonance  exists  among  the 
waves.  To  express  quantitively  the  nearness  of  the  above  resonances,  we 
introduce  the  two  detuning  parameters  ai  and  a2  defined  by 

u>3  - u>2  + ix>  1=  eoi , (36) 


Real(k3  - k2  + kj  = eo2. 
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where  ap  = 0(1).  Using  Eqs.  (36)  and  (37),  we  write 


- k2i)dx0. 


( (ill  + UJ3)To  = 02  + <<)  i 

(k!i  + 'ai 

(38) 

. + k . - k • )dx0 , 

2 ' 31  ll 

(39) 

i + k . - k ■ )dx0 , 

1'  2l  3T 

(40) 

where  kni-  stands  for  the  imaginary  part  of  kn  and 


<t>  = j* o2dx!-  oiT i (41 ) 

To  determine  the  An,  we  seek  a particular  solution  for  the  second- 
order  problem  in  the  form 

ui  = <hi(y;xi)exp(i0i)  + <h2(y;xi)exp(i02)  +ipi3(y;x!)  x 

exp(ie3)  + c.c. , (42) 

vi  = i^2i(y;xi)exp(i0i)  + ^22(y;xi)exp(i02)  + \Jj2 3 (y ; xa ) x 

exp(i03)  + c.c. , (43) 

Pi  = ^3i(y;xi)exp(i0i)  + i^32(y;xi )exp(i02)  + if>33(y;xi)  * 

exp(i03)  + c.c.  (44) 

Substituting  Eqs.  (27)  - (30)  and  (38)  - (44)  into  Eqs.  (22)  - (26)  and 
equating  the  coefficents  of  exp(i0i),  exp(i02),  and  exp(i03)  on  both 
sides,  we  obtain  the  following  equations: 

Mi  (*  k.)  = d .,  (45) 

lJ  2 J J lj 
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(46) 


M2(<l>  ,»  <P  ip  k.,  u.)  = d ., 
lj  2j  3j  J J 2 J 

M3 ip  ip  ii  kjt  W:)  = d s,  (47) 

lJ  ? J 3 J J J 3 J 

<P  i = <P  i = 0 at  y = 0,  (48) 

1 J 2 J 

ip  ^ > ip  j -*•  0 as  y -*■  °°«  (49) 

ij  2 j 

for  j = 1,  2,  and  3,  where  the  d. . are  given  in  Appendix  A. 

* J 

C.  Adjoint  Problem 

Since  the  homogeneous  parts  of  Eqs.  (45)  - (49)  are  the  same  as 
Eqs.  (31)  - (35)  and  since  the  latter  have  a nontrivial  solution,  the 
inhomogeneous  equations  (45)  - (49)  have  a solution  if,  and  only  if,  the 
inhomogeneous1  parts  are  orthogonal  to  every  solution  of  the  adjoint 
homogeneous  problem;  that  is. 


/ 


+ d,j!Tj  + d,j4tj)dy  = 0 for  J ■ ’■  2-  and  3l 


where  the  c*'s  are  the  solutions  of 

ki>  5 (Vo  - °Sj  = °- 

ko’  “j>  e i(u*ki  - “j’sj  * ,Ci0U» 

- 1 <°‘  - kKj  ■ °- 

s*j’  kr  “j>  5 1(u*kj  • “j’Cj  * 1kjsj 

1 


(50) 


(51) 


(52) 


‘ i (°2  * kj^2j  = °’ 


★ * 

^2 j " C3j 

= 0 

at 

0 

11 

>> 

* * 

C,J*  5,J 

- 0 

as 

y -*•  00 

(53) 

(54) 

(55) 
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- ■ — — 


Substituting  for  the  d. . from  Appendix  A into  Eq.  (50)  and  defining 

1 J 

x 

3j  = Ajexp[  f kj.dx],  (56) 

*x 

; we  obtain  the  following  differential  equations  for  the  evolution 

of  ai , a2 , and  a3: 

tr + Ir  = (ei  tt  - kn)ai  + e rr1'  »*«*«xp(-u).  (57) 

k %T  * ■ <*>  - k,,><»  ♦ E Tp  «.a,«xp(1*),  (58) 

l;  fr + w-  * (Ei  7f  - k,i)a>  * E (so 

where  u>^  = dwn/dkn  is  the  group  velocity  and  <j>  is  defined  in  Eq.  (41).  For 

spatial  modulation  only,  oi  = 0 and  9an/9t  = 0;  all  the  calculations 
presented  in  this  paper  are  for  this  case.  We  note  that  Eqs.  (57)  - 
(59)  account  for  the  combined  effects  of  the  nonparallelism  (i.e., 
growth  of  the  boundary  layer)  and  the  nonlinear  interaction.  If 
e « ei,  the  nonlinear  interactions  can  be  neglected  and  the  spatial 
variations  in  Eqs.  (57)  - (59)  reduce  to  the  nonparallel  solutions  of 
Refs.  11  and  12.  When  ei  « e,  the  effects  of  the  nonparallelism  are 
negligible;  that  is,  one  can  set  ei  = 0 and  all  the  coefficients  in  Eqs. 
(57)  - (59)  can  be  treated  as  constants. 

0 
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IV.  HARMONIC  RESONANCE 


The  interaction  between  a fundamental  Tollmien-Schiichting  wave  and 
its  second  harmonic  is  analyzed  using  a procedure  similar  to  that 
outlined  in  the  previous  section.  In  this  case,  instead  of  Eqs.  (57)  - 
(59)  we  obtain 


X dai  + 

oil  3t 


+ e 


a2aiexp(-i<j)). 


(60) 


X 3a 2 3a2 
ta)2  3t  3x 


where  $ is  defined  in  Eq. 


- k . )a2  + c 
2 * 

(41)  and 


a fexp(  i 4>) , 


(61) 


eo2  = Real (k2  - 2ki),  eoi=  w2  - 2ui,  (62) 

and  fi,  f2,  hu»  h22,  hJ2,  and  h2i  are  given  in  Appendices  B and  C. 

For  spatial  modulation  only,  aj  = 0 and  3an/3t  =0.  All  the 
calculations  presented  in  this  paper  are  for  this  case. 


% 
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V.  COMPUTATION  PROCEDURE 


Solutions  of  First-  and  Second-Order  Problems 


The  same  procedure  is  followed  in  solving  the  first-  and  second- 
order  problems  for  both  harmonic  and  combination  resonances.  Therefore, 
only  the  computation  methodology  for  the  solution  of  the  first-order 
problem  for  the  first  mode  is  outlined  here. 

Equations  (31)  - (33)  are  expressed  as  a system  of  first-order 
differential  equation  in  the  form 

(63) 


dz  - p7 

dx  " Gz’ 


where  z is  a 4 x 1 matrix  with  the  elements 


zi  = sn(y;xi),  z2  = Dcn(y;xi),  z3  =c2i(y;xi),  z4  = cst(y;xi), 

(64) 

and  6 is  a 4 x 4 matrix;  its  elements  are  given  in  Appendix  D. 


We  start  the  integration  of  Eqs.  (63)  at  y = yg,  where  yg  is  larger 


than  the  boundary-layer  thickness.  Hence,  U0  = 1,  DU0  = 0,  and 
D2U0  = 0 at  yg.  Then  the  matrix  G has  constant  coefficients  at  y = y£ 
and  Eqs.  (63)  have  solutions  of  the  form 


zi  = l ciiexp(Xiy)  for  i = 1,  2,  3,  and  4, 
j=i  J J 


(65) 


where  the  c. . are  constants,  the  X's  are  the  solutions  of 

' s) 


|G  - XI | = 0, 

and  I is  the  identity  matrix.  Equation  (66)  has  the  roots 


(66) 


1 >2 


= t kIf  X = ± [k2  + i ( k j - o)i )R]  ^2, 


3 »4 


(67) 
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Two  of  these  roots  have  positive  real  parts  that  make  the  solution  grow 
exponentially  as  y -*•  hence,  they  must  be  discarded  to  satisfy  conditions 
(35).  This  leaves  two  linearly  independent  solutions  that  dec»y  ex- 
ponentially with  y.  i 

The  eigenvalues  are  not  known  a priori  and  must  be  determined  along 

4 

with  the  eigenfunctions.  For  given  values  of  wi  and  R,  we  guess  a value 

for  ki  and  integrate  Eqs.  (63)  from  yg  to  y = 0.  If  the  quessed  value 

of  k3  does  not  staisfy  the  boundary  conditions  at  y = 0,  ki  is  incremented 

by  using  a Newton-Raphson  scheme  and  the  procedure  is  repeated  until  the 

boundary  conditions  are  satisfied  to  within  a specified  accuracy.  The 

integration  is  done  by  using  a computer  code  developed  by  Scott  and 
14 

Watts  . This  technique  orthonormal izes  the  solution  of  the  set  of 
equations  whenever  a loss  of  independence  is  detected. 

B.  Solution  of  Adjoint  Problem 

The  solution  procedure  is  exactly  the  same  as  that  for  the  first- 
order  problem.  The  coefficients  of  the  z matrix  are 

zi  » c*i(y;xi),  z2  = Dc*i(y;xi),  z3  = c!i(y;xi), 

★ 

z 4 = Cn(y;xi)  (68) 

and  the  adjoint  problem  has  the  same  eigenvalues  as  the  first-order 
problem. 

C.  Solvability  Conditions 

The  calculations  are  repeated  at  different  streamwise  locations  to 
evaluate  f.,  h..,  k.,  and  the  other  interaction  integrals  for  a given 

J J J J 

frequency  along  the  x-axis.  A fourth-order  fixed  step-size  Runge-Kutta 
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integration  scheme  is  used  to  solve  either  Ecs.  (57)  - (59)  for  com- 
bination resonances  or  Eqs.  (60)  and  (61)  for  harmonic  resonances  to 
find  the  amplitudes  of  the  waves  for  different  initial  amplitudes  of  the 
respective  modes. 
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VI.  Results  and  Discussion 


N 


The  analysis  presented  in  this  paper  is  applicable  to  both  two-and 
three-wave  interactions.  First,  we  present  and  discuss  numerical 
results  for  the  case  of  two-wave  interactions.  Then,  we  present  and 
discuss  numerical  results  for  the  interaction  of  three  waves  whose 
frequencies  are  such  that  F3  = F2  - Fa. 

A.  Two-Wave  Interactions 

The  numerical  results  presented  in  Ref.  6 show  that  the  amplitude 

of  a wave  a?  = ea2  must  exceed  a critical  value  before  it  can  generate 

and  amplify  its  subharmonic.  For  the  Blasius  flow,  the  critical  value 

is  approximately  29%  of  the  mean  flow.  This  is  for  the  case  when  the 

subharmonic  wave  has  an  infinitesimal  amplitude.  When  the  amplitude 
★ 

ai  = eai  of  the  subharmonic  wave  is  not  infinitesimal,  its  influence  on 
★ 

a2  should  be  taken  into  account.  The  equations  governing  this  influence 

are  Eqs.  (60)  and  (61)  whose  general  solution  is  not  available  yet.  The 

6 * 

previous  results  of  the  parametric  instability  model  show  that  ax 

★ 

oscillates  about  its  non-interaction  value  until  a2  reaches  the  critical 

value.  Figures  1 and  2,  obtained  by  numerically  solving  Eqs.  (60)  and 

* * 

(61),  agree  with  this  conclusion.  Initially,  a2  increases  while  ai 

oscillates  around  its  non-interaction  value. 

★ 

At  R ~ 580,  Fig.  1 shows  that  ai  starts  to  deviate  sharply  from  its 

★ 

non-interaction  value,  while  it  follows  from  Fig.  2 that  £na2  ~ -1.25 

★ ★ 
or  a2  = 0.286  at  this  location.  Hence,  when  a2  is  less  than  this 

★ 

critical  value,  aj  can  be  approximated  by  its  non-interaction  value; 
that  is 
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A 


I 


* * 

ai  = ajoexpt-k^.x  + ix)  (69) 

* 

where  ai0  and  i are  the  initial  amplitude  and  phase  of  the  subharmonic 

wave,  respectively.  If  we  substitute  Eq.  (69)  into  Eq.  (61)  and  neglect 

the  nonparallel  effects,  we  obtain 
* 

gj2-  + k^.at  = jj^-  at§exp[-(2kii  + ieo2)x  + 2ix].  (70) 

The  solution  of  Eq.  (70)  that  satisfies  the  initial  condition 
* * 

a2  = a2(0)  at  x = 0 can  be  written  as 


* 

a2 


= [a*(0)  + 


ai§exp(2ix)]exp(-k 


al§exp[-(2ki . + iea2)x  + 2ix]. 

(71) 


Equation  (71)  represents  an  approximation  to  a2  as  long  as  it  is  less 


than  the  critical  value  needed  to  trigger  the  parametric  instability  in 


the  subharmonic  wave. 

Next,  we  consider  the  generation  and  amplification  of  a second- 
harmonic  wave  by  a fundamental  Tollmien-Schlichting  wave.  We  consider 
the  following  three  cases:  (i)  fundamental  wave  is  stable  while  its 
second  harmonic  is  unstable,  (ii)  fundamental  wave  is  unstable  while  its 
second  harmonic  is  stable,  (iii)  both  fundamental  and  second- harmonic 
waves  are  unstable. 

When  the  fundamental  wave  is  initially  stable  while  its  second 
★ 

harmonic  is  unstable,  ai  decays  until  it  reaches  the  unstable  region  and 

then  it  increases  as  shown  in  Fig.  1.  For  Reynolds  Number  less  than 
* 

560,  ai  oscillates  around  its  non-interaction  value,  implying  a small 

"k 

initial  influence  of  its  second  harmonic  on  it.  Thus,  a2  can  be  approx- 
imated initially  by  Eq.  (71).  Figure  2 shows  that  the  values  obtained 
from  Eq.  (71)  are  in  good  agreement  with  those  obtained  by  numerically 
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integrating  Eqs.  (60)  and  (61)  for  R £ 560.  After  a short  initial 
distance,  the  second  term  on  the  right-hand  side  of  Eq.  (71)  decays. 

’ ft 

Then,  a2  can  be  approximated  by 


★ 

a2s 


(0)  + 


h2 1 

f2(2k~ii  + ieo2) 


*2 


aj  0exp(2ix)]exp(-k 


(72) 


★ 

as  long  as  a2  is  less  than  the  critical  value. 

Hence,  the  effect  of  the  fundamental  wave  on  its  second  harmonic  is 

* 

to  increase  its  initial  amplitude.  However,  as  a2  attains  large  values, 

★ ★ 
it  strongly  influences  ai  which  in  turn  strongly  influences  a2.  The 

result  is  an  accelerated  instability. 

For  the  case  when  the  fundamental  wave  is  initially  unstable  while 

its  second  harmonic  is  stable,  we  performed  calculations  for  waves  with 

the  frequencies  Fj  = 46.5  x 10"®  and  F2  = 93  X 10'®.  The  fundamental 

wave  is  in  the  unstable  region  at  R = 950  where  the  calculations  are 

started.  Thus,  its  unstable  downstream  of  R = 950.  Figure  3,  obtained 

★ 

by  numerically  integrating  Eqs.  (60)  and  (61),  shows  that  ai  hardly 

★ 

deviates  from  its  non-interaction  value.  On  the  other  hand,  a2  in- 
creases many  orders  of  magnitude  even  for  small  initial  amplitudes  of 
the  fundamental  wave  as  shown  in  Fig.  4.  In  these  calculations,  the 
initial  amplitude  of  the  second-harmonic  wave  is  taken  to  be  0.1%  while 
the  initial  amplitudes  of  the  fundamental  wave  are  0.1%  and  0.5%.  Since 

"ft 

hardly  deviates  from  its  non-interaction  value,  Eq.  (71)  is  expected 

★ 

to  be  a good  approximation  to  a2.  Figure  4 shows  that  the  values  obtained 

from  Eq.  (71)  oscillate  about  those  obtained  by  numerically  integrating 

Eqs.  (60)  and  (61).  Since  the  initial  values  are  very  small,  a2 

★ 

does  not  reach  the  critical  value  to  influence  ai.  After  a short  initial 
distance,  the  first  term  on  the  right-hand  side  of  Eq.  (71)  decays  and 
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a2  can  be  approximated  by 


yc  *■ 

ai 0exp[-(2k  ^ + ieo2)x  + 2ix], 


as  long  as  a2  is  less  than  the  critical  value.  Consequently,  the  effect 

of  the  interaction  is  to  produce  a second-harmonic  wave  that  grows 

approximately  at  a rate  that  is  twice  that  of  the  fundamental  wave. 

For  the  case  when  both  waves  are  unstable,  we  performed  numerical 

calculations  for  waves  having  the  frequencies  Fi  = 52  x 10"®  and 

F2  = 104  x 10"®  starting  near  R ~ 600.  Figure  5 shows  that  initially 

ai  deviates  slightly  from  its  non-interaction  values.  Hence,  Eq.  (71)  is 

★ 

expected  to  be  initially  a good  approximation  to  a2.  Figure  6 shows 

that  the  numerical  values  obtained  from  Eq.  (71)  are  in  good  agreement 

with  those  obtained  by  numerically  integrating  Eqs.  (60)  and  (61)  when 
★ 

a2  is  less  than  0.29.  Thus,  in  this  case,  the  effect  of  the  interaction 
on  the  second-harmonic  wave  is  to  increase  its  initial  amplitude  and  to 
produce  a term  that  grows  at  a rate  that  is  twice  the  growth  rate  of 
the  fundamental  wave.  Due  to  the  fact  that  both  waves  are  initially 
unstable,  the  interaction  is  more  effective  in  this  case  than  in  the 
preceding  two  cases. 

B.  Three-Wave  Interactions 

9 8 10 

In  their  experimental  studies,  Kachanov,  et  al  , Miksad  , Norman  , 

and  Sato^  introduced  two  separate  waves  of  different  frequencies  into 

the  flow  that  was  being  studied.  They  observed  the  growth  of  a wave 

whose  frequency  is  equal  to  the  difference  frequency.  Kachanov  et  al 

-fi  — fi 

used  the  frequency  pairs  Fi  = 88  x 10  and  F2  = 104  x 10  and 
-fi  -fi 

Fa  = 88  x 10  and  F2  = 120  x 10"  to  analyze  the  growth  of  the  associated 
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-6 

difference-harmonic  waves  at  F3  = 16  x 10  (i.e.  F3  = F2  - F3)  and 

F3  = 32  x 10"®  in  a boundary-layer  flow  over  a flat  plate.  Using  the 
same  frequency  pairs,  we  determined  the  amplitudes  of  the  fundamental 

k k k 

waves  (ai  and  a2)  and  the  difference-harmonic  wave  a3  by  numerically 
solving  Eqs.  (57)  - (59).  Figures  7 and  8 show  the  large  increases  in 
the  amplitudes  of  the  above  difference-harmonic  waves  due  to  the  nonlinear 
interaction  of  the  waves  with  the  frequencies  Fj  and  F2.  The  amplitude 
of  the  difference-harmonic  wave  at  R ~ 715  is  increased  120-200  times 
its  non-interaction  value. 

However,  the  amplitudes  of  the  fundamental  waves  change  very  little 
from  their  non-interaction  values  as  shown  in  Fig.  9 for  Fj  = 88  x 10"®. 

★ k 

Thus,  ai  and  a2  can  be  approximated  by 

k k k k 

aj  = a3 0exp(-k  .x  + in),  a2  = a20exp(-k  -x  + ix2),  (74) 

1 1 2 1 

* * 

where  ai0  and  a20  are  the  initial  amplitudes  and  n and  t2  are  the 
initial  phases  of  the  fundamental  waves.  Substituting  Eq.  (74)  into 
Eq.  (59)  and  neglecting  the  nonparallel  effects,  we  obtain 


+ k 4 at 

31 


★ * 
a 1 oa2  0 


exp[-(ki . 


+ k . + iea2)x 
21 


+ i(xi  + t2)]  (75) 

The  solution  of  Eq.  (75)  that  satisfies  the  initial  condition 

★ k 

a3  = a3(0)  at  x = 0 can  be  expressed  as 

at  = [at  (0)  + •f3-(k-/'|1t  7 +~~i ea2 ) atoatoexp(in  + ix2 )]  x 

ii  21 


exp(-kMx)  - m . u-  ) a!„atoexp[-(k  , a k , 

1 ' 2 1 

+ iea2)x  + i (x  j + x2)].  (76) 

k 

Equation  (76)  represents  an  approximation  to  a3  as  long  as  the  amplitudes 
of  the  fundamental  waves  do  not  deviate  from  their  non-interaction  values. 
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Figures  7 and  8 show  that  Eq.  (76}  is  initially  in  good  agreement  with 
those  obtained  by  numerically  integrating  Eqs.  (57)  - (59). 

According  to  Eq.  (76),  the  difference-harmonic  wave  grows  at  a rate 
that  is  the  sum  of  the  growth  rates  of  the  fundamental  waves.  Since  the 
fundamental  waves  are  unstable  at  R = 430,  where  the  calculations  are 
started,  the  difference-harmonic  wave  amplifies  considerably,  inspite  of 
the  fact  that  it  is  stable  in  the  absence  of  the  interaction.  These 
results  are  in  qualitative  agreement  with  the  experimental  observations 
of  Refs.  7-10. 
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APPENDIX  A 


3A.  3CH 

dlJ  = ( 3*7  Clj  + Aj  37T  K 


d2j  = d2j0  + d2ji , 

o 3A.:  3^i. 

d2j0=  "(UoCi j + S3j  - -r  kj^j  - VoD^j)  a*f  + c - u°  a*T 


+ k1  (k 


3Cn  dk. 


3xi  R ' j 3Xj  dxi  j 3xj  ^ 1 j J j 


CiJA,, 


d2u  = [i(k2  - ka)Ci2^i3  + C22DCi3  + D£i2£23]A2A3  x 

exp[-i<j)  - J( k2.  + k3i  - k!.)dx], 

d22i  = [i (kj  + k3)ci iCi 3 + €2iD£i3  + D^j iC23]AiA3  x 


exp[  i 4> 


H - J (kli  + k31  - k2.)dx]. 


d23i  = (“ikjCiiCi2  + ik2Ci2^x i + C22D^n  + D^i2C2i)A2^i  x 


expC-i  <J>  -J (kXl  + k2l  - k3.)dx]. 


d3j  = d3j°  + d3j*’ 


04  0*4  0^2  -4 

d3j0=  - (UoC2j  - r kjC2j)  -fii  - [Uo  3^  + V0DC2j  + c* jDVo 


94  dk.  3^2 4 

• R-  ( d^  + kj  3771  )]Aj’ 


= (-ik3Cl2C2  3 + ik2C22Cl 3 + C22DC23  + 2 2 £ 2 3 )A2A3  X 


exp[-i<f,  -/(k2i  + k3i  - kli)dx]. 


= (^k3^i 1C23  + ikj^21^13  + C2lD^23  + D^21^23)A1A3  X 


exp[id>  - f( kn  + k3.  - k2 . )dx] , 


(A10) 


APPENDIX  B 


? { * 2ik.-  * 2ik . 

J0  ( " C(U°  ’ * )Clj  + C3j]C2J  ’ (U°  ' ^ 


OjttjJdy 

if?  4 


( j - 1,  2,  3), 


3Cii  3^3 

+ cu' 


3^» 

. - — - ( k * 

J R v j 3Xi 


ur4  j\i»  if  2ik . 3^2^ 

m Clj}  + 377  Ctjfcaj  + C(Uo  - V-  } 37^-  +VoDC2 


3Xi  'u“^j 


o-j  * ) 

*jDV°  ■ T dxf  ^2j]c3j  j dy  U = ]*  2»  3>» 


■7! 


[i (k2  ■ k3)Cl2?13  + C22DC13  + D^ 1 2 £ 2 3 2 1 


+ ( -ik3Sl2C23  + ik2C22^ 13  + C22DC23  + D?22C23)C 


■7!1' 


tij  dy. 


[i(ki  + k3)Cjj^i3  + C21DC13  + DCuC23]?22 


+ (ik3CllC23  + iklC2lSl3  + C2lD^23  + 0^21^23)^32 


* ) 

32>dy, 


■/!” 


[i(k2  - klklzCll  + C22DCu  + DCl2^2l]C23 


+ (-ikl^2lCl2  + i k2^22Cl 1 + C22DC2I  + D^22  C2 1 ) £3 


3 3 j dy . 
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APPENDIX  C 


'12 


oo 


u + CaiD^iz  + i(k2  - ki)c12rjU]c*i 


+ [1k2C22Cn  - 11ciC21Ci2+  c220c2,  + C21Dc22]cti|  dy  (Cl) 

00 

*>n  -j(j[(ikiCii  ♦iiiDtnkt*  ♦ (IkiCnCi,  + C!iDC!i)it2])dy  (C2) 


])C 


APPENDIX  D 


9n  = 0 > 9w  * 1»  9i3  = 0*  9i4  = 0 (D1) 

921  = i(U0k  - w)R  + k2,  922  = 0,  g23  = R 924  - ikR  (D2) 
9si  55  - ik,  g32  = g33  = g34  = 0 (D3) 

9m  ■ 0,  g„2  = -ik/R,  g43  = -[i(U0k  - w)  + k2/R],  g44  = 0 (D4) 

Q is  a 2 x 4 matrix  consisting  of  the  last  two  rows  of  the  matrix 


B"1.  The  matrix  B has  the  elements: 

bi i = bi2  = bj3  = bi4  = 1 (D5) 

b23  = -k,  b22  = k,  b23  = k,  b24  = - k (D6) 

b3i  = i,  b32  = ik/k,  b33  = -i,  b34  = -ik/k  (D7) 

b4i  = (co/ k - 1),  b42  = 0,  b43  = (u/k  - 1),  b44  = 0 (D8) 

where 

k = [k2  + i(k  - o))R]1/z  (D9) 


subharmonic  resonance  with  a fundamental  wave 


involved  In  resonance  with  its  subharmonic  at  F 


nvolved  In  resonance  with  its  harmonic  at  F 


when  involved  in  resonance  with  a wave  at  F.  = 46.5  x 10" 


Figure  8.  Amplitude  of  the  difference-harmonic  wave  a 
generated  from  a combination  resonance  of  F 
and  F2  = 104  x 10‘6. 


